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Abstract —Air pollution monitoring is a very popular research 
topic and many monitoring systems have been developed. In this 
paper, we formulate the Bus Sensor Deployment Problem (BSDP) 
to select the bus routes on which sensors are deployed, and 
we use Chemical Reaction Optimization (CRO) to solve BSDP. 
CRO is a recently proposed metaheuristic designed to solve a 
wide range of optimization problems. Using the real world data, 
namely Hong Kong Island bus route data, we perform a series of 
simulations and the results show that CRO is capable of solving 
this optimization problem efficiently. 

Index Terms —Air pollution, public transportation, evolution¬ 
ary algorithm, chemical reaction optimization. 

I. Introduction 

A ir pollution has raised great concern over the past few 
decades due to the increasing expansion of industries. It 
has caused many serious problems, including climate change, 
loss of biodiversity, changes in hydrological systems, acid rain, 
and stress on the system of food production (TJ. It is also 
known that some of the chemical pollutants in the air can 
increase the occurrence of diseases such as lung cancer and 
pneumonia @0. Since hazardous gases can spread over a 
very large region and causing huge and irreparable damage 
a, there is a growing demand for air pollution monitoring 
systems. 

Many air pollution monitoring systems utilizing smart 
sensor networks and wireless systems have been proposed 
in the recent literature 00. But most of these systems 
use individually designed facilities to collect and transmit 
detected data, and the facilities are installed on stationary 
bases. These systems require a large number of sensors in 
order to provide a satisfactory coverage of the whole area, 
rendering them very expensive. To increase the coverage with 
a limited number of sensors, one can have embedded sensors 
installed on moving objects, e.g., a vehicle or an animal. With 
the movement of these objects, a large area can be covered. 
Since the air pollution condition changes relatively slow and 
can be regarded as constant in a short period of time, the 
mobile sensor system can achieve a larger coverage without 
losing too much accuracy. Most previous work focuses on 
implementation aspects, e.g. the design of mobile sensor and 
installation on the bus 00. However, no literature has been 


published concerning the selection of bus routes on which 
sensors are to be deployed. In this paper, we propose a novel 
optimization problem solving this selection problem. 

To select the bus routes on which sensors are to be deployed 
so as to minimize the number of sensors required for a 
satisfactory coverage, we formulate an optimization problem 
called Bus Sensor Deployment Problem (BSDP). To solve 
this problem, we employ the Chemical Reaction Optimization 
(CRO) technique. CRO is a population-based general-purpose 
optimization metaheuristic which mimics the transition and 
interaction of molecules in a chemical reaction. In chemical 
reaction there is a natural tendency for the potential energy 
of the reactant molecules to decrease until it reaches a stable 
energy state 0- CRO utilizes this tendency to guide molecules 
to explore the solution space and to find the global minimum. 

In our model, we assume that the whole area is divided 
into square grids of the same size and we assume that the 
air pollution condition in the same grid are similar. The route 
of a bus are divided into segments according to the boundary 
of the grids and we consider the mid-point of these segments 
as the sensing point. When the bus completes its route, the 
stored data are uploaded wirelessly to the base station at the 
bus terminus. 

The rest of this paper is organized as follows. Section II 
reviews the related work on air pollution monitoring system. 
The problem to be solved using CRO is described in Section 
III, followed by a detailed framework and algorithm design in 
Section IV. The simulation results are reported and discussed 
in Section V. Finally we conclude the paper in Section VI. 

II. Related Work 

Air pollution monitoring is a hot research topic due to 
the increasing concern on the adverse effects of pollution. 
Kularatna et al. proposed an environmental air pollution 
monitoring system in 0 focusing on CO, N0 2 , and S0 2 
detection. The proposed system is based on a smart sensor 
converter installed with an application processor which can 
download the pollution condition for further processing. Tsow 
et al. proposed a wearable and wireless sensor system for 
real-time monitoring of toxic environmental volatile organic 


compounds in m. Jung et al. proposed an air pollution geo¬ 
sensor network to monitor several air pollutants in 0. The 
system consists of 24 sensors and 10 routers, and provides 
alarm messages depending on the detected pollutants. Gao 
et al. proposed a wireless mesh network to cover a given 
geographic area using embedded microprocessors consisting 
of sensors and wireless communication in 0- Kwon et al. 
proposed another outdoor air pollution monitoring system in 
eh This system uses ZigBee networks to transmit the sensed 
pollutant density levels. The above systems are all air pollution 
systems utilizing mobile sensors to achieve high coverage, 
but they all need proprietary equipment to accommodate the 
movement requirements of the system. Gil-Castineira et al. 
proposed an air pollution detection system based on the public 
transportation system and tested it in a small scale experiment 
0. However, there is no sensor deployment algorithm that 
can efficiently utilize the available resources. So in this paper, 
we formulate an optimization problem to deploy the sensors 
so as to utilize them efficiently. The optimization problem is 
solved with CRO. 

Many optimization problems have been solved using CRO 
since □2). Xu et al. used CRO to solve task scheduling 
problem in grid computing ini. This problem is a multi¬ 
objective NP-hard optimization problem. Lam et al. proposed 
a population transition problem in P2P live streaming and 
solved this problem using CRO in d. Lam and Li also 
solved the cognitive radio spectrum allocation problem in 
m. Several variants of CRO were proposed to solve the 
optimization problem and a self-adaptive scheme was used 
to control the convergence speed of CRO eh Yu et al. pro¬ 
posed a CROANN algorithm from real-coded version of CRO 
ESI to train artificial neural networks (ANNs). CROANN 
used a novel stopping criteria to prevent the ANNs from 
being over-trained and the simulation results demonstrated that 
CROANN outperformed most previously proposed EA-based 
ANNs training methods as well as some sophisticated heuristic 
training methods. This shows that CRO has great potential to 
tackle different optimization problems like BSDP discussed in 
this paper. 


III. Problem Formulation 

In this paper, we formulate BSDP as a grid coverage 
problem. The monitoring area A is divided into p x q grids. 
There are n bus routes and each bus route li,. i = 1,2,..., n 
passes through multiple grids. If a bus route passes through 
a grid, then we say this bus route “covers” the grid. We 
then define a coverage threshold c such that a grid is “fully- 
covered” if it can be covered by c or more routes. A solution S 
is given by a vector S = [si, s 2 ,..., s n ], where s, = 1 means 
there is a sensor installed on the buses for Ri and ,s, = 0 
means no sensors installed for It, . With all the bus routes’ 
sensor-installation information .s t collected, we can make a 
coverage graph Gs f where the solution Sf = [1,1,. ..,1] 
means all bus routes are equipped with sensors, and we define 
the total number of covered grids in Gs f is ts f ■ Then given 
a random solution S' we can compute its coverage graph Gs> 



and its covered grid number ts' . It is natural to assume that the 
pollutant level in each grid remains constant in a short period 
of time like one hour 0 and we set the sensing interval for 
each grid to be one hour. Since all buses we use for simulation 
can finish their routes in this period of time, we consider that 
there is only one bus with sensor installed running on the 
same route simultaneously, and the bus can pass all the grids 
in each sensing interval. Take Fig. 1 as an example. The two 
black circles are two bus routes. The 28 grids with shading 
are covered by at least one route. The four grids in the center 
with darker shading are covered by two routes. So in this 
example if we set c = 1 or 2, ts f = 28 or 4, respectively. In 
later calculation we only consider these tg f grids. If we only 
install sensors on buses running on the left circle, the ts' is 
16 or 4 with c = 1 or 2. If we desire to have a shorter sensing 
interval, say to times instead of once per hour, we can deploy 
sensors on m of the buses on each route, with bus start times 
of 1 /to hours apart. 

In this problem we have two major factors to evaluate a solu¬ 
tion: the coverage percentage and the total sensor number. The 
former concerns the performance as well as data accuracy of 
the system while the latter concerns the expense of deploying 
the system. Since in different scenarios the system may have 
different requirements, we introduce a weight coefficient a to 
balance the two factors. The objective function for BSDP is 
given as follows: 

min (1-— ) x a + ^~' z=1 S ‘ x (1 — a). (1) 

t Sf n 

This objective function is composed of two parts: the 
percentage of uncovered grids over the total coverable grids 

t, y n _ si 

and the percentage of sensor deployment this 

paper, we adopt a = 0.5 in our simulations in Section V to 
simulate a general case that both the coverage and the cost of 
purchasing sensors are important. 





























IV. Algorithm Design 

In this section, the detailed design of our algorithm to solve 
BSDP is given. First we briefly introduce how CRO works, 
and then the encoding scheme and operators employed for 
optimization are presented. 

A. Chemical Reaction Optimization 

CRO is a kind of variable-population-size-based meta¬ 
heuristics that can solve optimization problems efficiently. 
CRO exploits the natural tendency of chemical reactions to 
reduce the total potential energy in reactant molecules to 
search the solution space and to find the global optimum. 
In a chemical reaction process, the molecules with potential 
energy (PE) as well as kinetic energy (KE) are put into a 
closed container. When a collision happens, the molecules will 
change their structure to transform PE to KE or vice versa, or 
or just release the energy to the environment. If we consider 
the energy states of molecules as a surface, this procedure can 
be considered as molecules rolling down the energy surface 
to reach the lowest energy state. CRO utilizes this natural 
tendency to perform optimization. 

In CRO, each molecule has a molecular structure u> and 
two kinds of energy, i.e. KE and PE. The molecular structure 
stands for a feasible solution to the problem, PE is the 
objective function value for the solution, and KE is set as 
a tolerance for the molecule to move to another energy state 
with higher energy. We use four different types of elementary 
reactions to imitate all kinds of molecular collisions, namely, 
on-wall ineffective collision, decomposition reaction, inter- 
molecular ineffective collision, and synthesis reaction. The 
four elementary reactions cooperate with each other to search 
the minimums while maintaining a wide population diversity. 

When CRO algorithm starts, some randomly generated 
molecules are initially put into a closed container. Then in 
each iteration one collision takes place in the container. The 
collision can either be one molecule colliding on the wall, 
or two molecules colliding with each other. We divide the 
four different elementary reactions according to the molecules 
involved in the collision into two categories: uni-molecular 
collisions and inter-molecular collisions, and we first randomly 
select one category. The former includes the on-wall ineffec¬ 
tive collision and decomposition, while the latter includes the 
inter-molecular ineffective collision and synthesis. After the 
reaction category has been decided, the system will randomly 
select molecule(s) to participate in the chemical reaction. The 
system then check the energy of the selected molecule(s) to 
determine which exact elementary reaction shall happen and 
then the corresponding operator is performed on the involved 
molecule(s). The final step before the end of each iteration 
is the performance check. The objective function value(s) of 
the newly generated molecule(s) is computed and compared 
with previous value(s). If the new value(s) can satisfy the 
energy conservation conditions discussed in G2. the new 
molecule(s) is accepted and substitutes the original molecules. 
Otherwise the new molecule(s) is discarded. This completes 
an iteration of CRO. After the number of iterations reaches a 


certain number or other stopping criteria is met, the algorithm 
terminates. Interested readers can refer to ED and H6) for 
detailed description of the algorithm as well as its pseudocode. 

B. Encoding Scheme 

We use a simple encoding scheme to encode a solution 
of BSDR A solution is formulated as a vector of n binary 
numbers. Each element in the solution stands for whether to 
install sensors for the specific route or not. 

C. Operators 

CRO has four different types of elementary reactions, which 
correspond to different functionalities. So we design a corre¬ 
sponding operator for each of them. We also design an initial 
solution generator to generate the solution structures of new 
molecules. We generate all random numbers uniformly in the 
given ranges, unless stated otherwise. 

1) Initial Solution Generator: This initial solution genera¬ 
tor is designed to generate new molecular structures, which is 
triggered when CRO is initialized or a decomposition happens 
oa. We randomly assign 0 or 1 to each element in the 
vector to generate new molecules. Its pseudocode is given in 
Algorithm 1 below: 


Algorithm 1 InitialGen (w) 

1: for all Elements e in u> do 

2: Randomly generate a real number n £ [0,1). 

3: if n < 0.5 then 

4: e = 0 

5: else 

6 : e = 1 

7: end if 

B: end for 


2) Neighborhood Search Operator: This operator is applied 
to the two ineffective reactions, namely the on-wall ineffective 
collision and the inter-molecular ineffective collision. It is 
designed to generate a new molecular structure oj' from the 
neighborhood of the given molecular structure to. Its main 
purpose is to perform a detailed local search for potentially 
better solutions E2. A random-toggle scheme is used to 
perform this operation. We first randomly pick an element 
c, from ui and then update the value by e,; = 1 - e, . This 
operation can efficiently perform a neighborhood search on 
the solution space without losing accuracy. The pseudocode 
of this operator is given in Algorithm 2 below: 

Algorithm 2 INEFFECTIVE (w) 

1: Generate a random integer i smaller than the total number 
of elements in a solution 
2: Find the i th element in ui 

3: — 1 6i 


In an on-wall ineffective collisions, one molecule is in¬ 
volved and this operator can be directly applied to change the 
molecular structure. However, in an inter-molecular ineffective 











collision, since two molecules are involved, we separately 
manipulate them using this operator to perform local search. 

3) Decomposition: This operator is used for generating two 
new molecular structures wj and uj' 2 from the given molecular 
structure u>. This operator mainly focuses on helping the 
algorithm to jump out of local minimums na by making 
severe changes with energy sharing. The pseudocode of this 
operator is given in Algorithm 3 on the next page. 


Algorithm 3 Decomposition (w) 

1: Copy ui to wj and u>' 2 

2: for all Elements e in wj’s and ui' 2 ’ s molecular structure 

do 

3: Randomly generate a real number n £ [0,1). 

4: if n < 0.5 then 

5: e = 1 — e 

6: end if 

7: end for 


The original molecule is copied to two new molecules and 
each element in the new molecular structure is individually 
toggled with a probability of 0.5. If this operation can satisfy 
the energy balance rule mentioned in lfl2l then this reaction is 
accepted, i.e. the original molecule is discarded and the two 
new molecules are put into the container. 

4) Synthesis: This operator is used for generating one new 
molecular structure u>' from two given molecular structure 
u >i and uj 2 . This operator can perform a general local search 
while preventing the molecules from being stuck in the local 
minimums sm The pseudocode of this operator is given in 
Algorithm 4 below: 


Algorithm 4 SYNTHESIS (wi, w 2 ) 

1: for all Elements e in u/’s molecular structure do 
2: Randomly generate a real number n £ [0,1). 

3: if r > 0.5 then 

4: e counterpart in u>\ 

5: else 

6: e counterpart in oj 2 

7: end if 

8: end for 


In this operator, the new molecule is composed of the two 
given original molecules and each element of the solution 
is equally likely to be selected from each of the original 
molecules. 

V. Simulation Results 

In this section we will first introduce the data used for 
simulation of BSDR Then the detailed simulation parameter 
settings, results, and comparisons are presented. 

A. Simulation Data 

In order to make the simulation results persuasive, we 
utilize the real data of the Hong Kong Island bus routes in 
the simulation. We selected 91 bus routes from the Citybus 



Fig. 2. Grids and 91 Selected Routes for Simulation and 


transportation system 03 which covers most of the accessible 
areas on Hong Kong Island. There are totally 2277 stations 
for the 91 bus routes and for simplicity, but without losing 
generality, adjacent stations on the same routes are directly 
connected and buses are supposed to run on the connected 
graph. The plot of all routes is given in Fig. 2. 

The whole area is divided into 16 x 11 = 176 1 km x 1 km 
grids. For different coverage threshold c we have different 
ts r For instance, if c = 5 then ts f = 40. The reason why 
ts f cannot reach 176 is that there are no bus routes passing 
through some rural grids. 

B. Analysis of CRO parameters selection 

The performance of CRO is greatly influenced by the proper 
selection of the optimization parameters HI. The ratio of 
occurrence of different reactions, the tolerance of molecules to 
jump to a high energy state and the energy consumption rate 
are all key factors that drive significant impact on the quality 
of the final solutions. So it is essential to analyze and select 
a proper combination of parameters for the simulation. In this 
analysis, we use the Hong Kong Island bus route data and 
set c = 5. The results of testing the parameters are generated 
by computing the average of 50 trials with 10 000 function 
evaluations. Here are brief introductions to the six parameters 
and interested reader can refer to 02 and na for more 
information. 

1) EnBujf: This parameter describes the initial energy 
buffer size of the container. When the algorithm iterates, this 
energy can be transfered from and to the molecules’ PE or 
KE. 

2) IniKE: This parameter describes the initial kinetic en¬ 
ergy each molecule holds. It decides the tolerance of accepting 
bad molecules. 

3) CollRate: This parameter describes the fraction of an 
elementary reactions being an inter-molecular collision. It 
functions when an iteration starts and decide the category of 









































TABLE I 

Analysis on Performance for Different CRO Parameters 


EnBuff 

IniKE 

CollRate 

Value 

Mean 

Std. 

Best 

Value 

Mean 

Std. 

Best 

Value 

Mean 

Std. 

Best 

0 

0.423302 

0.0863343 

0.398626 

0 

0.428764 

0.11517 

0.387637 

0.1 

0.429088 

0.0929551 

0.398626 

10 

0.422423 

0.102802 

0.395604 

10 

0.430577 

0.158681 

0.387637 

0.2 

0.429736 

0.111531 

0.387637 

100 

0.423159 

0.0828353 

0.395604 

100 

0.425835 

0.111766 

0.395604 

0.4 

0.428396 

0.0933636 

0.401648 

1000 

0.428346 

0.0879891 

0.398626 

1000 

0.427692 

0.112578 

0.398626 

0.6 

0.43139 

0.0969814 

0.409615 

5000 

0.428121 

0.0823935 

0.40467 

5000 

0.427286 

0.102947 

0.387637 

0.8 

0.426929 

0.0986975 

0.398626 

10000 

0.423099 

0.0932018 

0.37967 

10000 

0.431253 

0.108105 

0.398626 

0.9 

0.430945 

0.108942 

0.398626 

LossRate 

DecThres 

SynThres 

Value 

Mean 

Std. 

Best 

Value 

Mean 

Std. 

Best 

Value 

Mean 

Std. 

Best 

0.1 

0.432907 

0.125838 

0.398626 

100 

0.450401 

0.147143 

0.409615 

10 

0.411236 

0.0971248 

0.376648 

0.2 

0.43133 

0.127488 

0.384615 

300 

0.409709 

0.091254 

0.37967 

50 

0.413066 

0.0835925 

0.384615 

0.4 

0.423527 

0.109854 

0.387637 

500 

0.413154 

0.0904357 

0.387637 

100 

0.406973 

0.0803982 

0.387637 

0.6 

0.414357 

0.0957439 

0.384615 

1000 

0.410253 

0.0943316 

0.376648 

300 

0.410703 

0.094436 

0.376648 

0.8 

0.410626 

0.0908277 

0.376648 

3000 

0.419896 

0.121163 

0.37967 

500 

0.41239 

0.0972623 

0.387637 

0.9 

0.422995 

0.093694 

0.387637 

5000 

0.427819 

0.210808 

0.376648 

1000 

0.411368 

0.0983823 

0.387637 


TABLE II 
CRO Parameters 


Parameters 

Values 

Function Evaluation Limit 

10 000 

Initial Population Size 

20 

Initial Energy Buffer Size 

5000 

Initial Molecular Kinetic Energy 

5000 

Molecular Collision Rate 

0.4 

Kinetic Energy Loss Rate 

0.8 

Decomposition Threshold 

300 

Synthesis Threshold 

100 


reactions for this iteration. The larger this parameter is, the 
larger possibility that inter-molecular collision will happen. 

4) LossRate: This parameter describes the energy loss rate 
when an on-wall ineffective collision happens. 

5) DecThres: This parameter describes a threshold for 
decomposition reaction. This parameter was named “alpha” 
in m. When the iteration is decided to be a uni-molecular 
collision, the algorithm will decide whether the decomposition 
reaction shall be conducted or not by using this parameter. 

6) SynThres: Similar to DecThres, this parameter describes 
the happening threshold for synthesis reaction. This parameter 
was named “beta” in 02 and operates when inter-molecular 
collision is selected. 

From the analysis results listed in Table I we can see that 
although the performance of the simulation can stay at a 
relatively high level, proper parameters can improve the final 
performance in a small scale. So we select a combination of 
parameters whose performance in mean result as well as the 
standard derivation is relative good. The proper combination 
for simulation is listed in Table II. 

C. Analysis on the Impact of c to Result 

To analyze the impact of c on the BSDP optimization result, 
we adopt different values of c and perform simulations using 


the previously stated parameter combination. The results are 
listed in Table III and are generated by computing the average 
of 50 trials with 10 000 function evaluations. 

In Table III, the first column “Objective Function Result” 
presents the raw data obtained from the objective function 
listed in Section III. The “Full-coverage Result” column 
presents ts f described in Section III, with respect to different 
value of c. If the covered grid count is high, then the largest 
possible cover percentage over the whole area is also high. 
Since in the Hong Kong Island data there are many grids 
uncovered (because of no inhabitation or sea area), the per¬ 
centages are relatively low. The “Average Result” presents the 
mature data of the average value of the 50 best result generated 
from the trials. The “Best Result” presents the mature data for 
the best-performing result from the 50 best result. 

The selection of c has a large impact on the final result. 
Since c defines the coverage threshold, a larger c will result in 
more buses passing through a grid. However the total number 
of bus routes is a fixed number, and in some remote area it 
is possible that there are totally less than c bus routes. On the 
contrary, a smaller c can result in a possibly higher coverage 
rate, the collected data accuracy is less. So the selection of c 
shall be decided by the real-world requirement of the system. 

D. Comparison among Chemical Reaction Optimization, A 
Random Method, and Simple Genetic Algorithm 

We compare the performance of CRO with a greedy ap¬ 
proach called Simple Random Method (SRM), and with Sim¬ 
ple Genetic Algorithm (SGA). SRM generates 10 000 random 
solutions and stores the best-so-far solution. For each element 
in the solution SRM will perform the randomization elaborated 
below. First we generate a random number n £ [0,1) as 
“base”, then we generate another random number m £ [0,1) 
independently as “target”. If to > n then this element is set 
to 1; otherwise, the element is set to 0. The reason why we 
generate two random numbers instead of one for randomiza- 

























TABLE III 

Analysis of the Impact of c on the Results 



Objective Function Result 

Full-coverage 

Result 

Average Result 

Best Result 


Mean 

Std. 

Best 

Covered Grids 

Percentage 

Coverage 

Routes Count 

Coverage 

Routes Count 

1 

0.153231 

0.0291638 

0.142857 

65 

36.9318% 

99.7538% 

13.72 

100.000% 

13 

2 

0.215108 

0.0434328 

0.201872 

54 

30.6818% 

97.8889% 

17.70 

96.2963% 

15 

3 

0.278838 

0.0556753 

0.262166 

49 

27.8409% 

97.9184% 

23.48 

97.9592% 

22 

4 

0.359091 

0.0953412 

0.329670 

44 

25.0000% 

98.0909% 

30.94 

100.000% 

30 

5 

0.410451 

0.0836293 

0.387637 

40 

22.7273% 

95.9000% 

33.62 

97.5000% 

33 

6 

0.482267 

0.1541760 

0.452574 

38 

21.5909% 

91.1579% 

35.84 

92.1053% 

34 


TABLE IV 

Comparison among CRO, SRM, and SGA 



CRO Result 

SRM Result 

SGA Result 


Mean 

Std. 

Best 

Mean 

Std. 

Best 

Mean 

Std. 

Best 

1 

0.153231 

0.0291638 

0.142857 

0.401363 

0.119488 

0.389011 

0.40633 

0.150758 

0.323077 

2 

0.215108 

0.0434328 

0.201872 

0.450981 

0.113744 

0.455230 

0.455551 

0.136338 

0.407204 

3 

0.278838 

0.0556753 

0.262166 

0.508006 

0.104393 

0.475667 

0.522512 

0.127871 

0.459969 

4 

0.359091 

0.0953412 

0.329670 

0.55979 

0.095316 

0.563437 

0.567597 

0.084297 

0.53047 

5 

0.410451 

0.0836293 

0.387637 

0.587005 

0.110423 

0.575549 

0.599115 

0.107698 

0.561538 

6 

0.482267 

0.1541760 

0.452574 

0.668028 

0.129847 

0.663389 

0.670804 

0.132876 

0.621747 


tion is that we cannot set a fixed threshold controlling the ratio 
between sensor-equipped routes and other routes. But with this 
double-randomization technique, the elements are individually 
generated and all have different probabilities of being 1. The 
pseudocode of SRM is given in Algorithm 5 below: 


Algorithm 5 SRM 

1: Set GlobalMin to be a large number 

2: while Function evaluation count is not exhausted do 

3: Initiate a new solution s 

4: for all Elements e in s do 

5: Randomly generate a real number n£ [0,1). 

6: Randomly generate another real number m £ 

[ 0 , 1 ). 

if m > n then 

e = 1 

else 

e = 0 

end if 
end for 

if The objective function value of s is smaller than 
GlobalMin then 
14: GlobalMin = objective function value of s 

15: end if 

16: end while 


B 

9 

10 

11 

12 

13 


We follow |(T9l to program SGA and we adopt the crossover 
rate as 0.5 and permutation rate as 0.1 for this simulation. For 
all three algorithms, we set the function evaluation limits to 10 
000. Since the parameter c shall be decided by the real-world 
requirement, we compare the performance of CRO, SRM and 
SGA with c £ {1,2,3,4, 5, 6}, respectively. All the simulation 
results are presented in Table IV. 


From Table IV, we can see that CRO can outperform SRM 
and SGA dramatically with different values of c for the mean, 
the standard deviation and the best result. This shows that CRO 
is an efficient way in solving BDSP optimization problem. 

VI. Conclusion and Future Work 

In this paper, we propose a novel air pollution monitoring 
system by deploying sensors in a public transportation system. 
We also formulate a new optimization problem for selecting 
the buses to deploy the sensors, called BSDR The main idea 
is to install sensors on buses and with the movement of the 
buses, the sensors can cover a much larger area compared 
with stationary sensor stations. This raises the problem of 
selecting buses to install sensors. We use CRO to solve 
this optimization problem since CRO has been shown to 
be powerful in optimizing similar problems mm. In our 
simulation, we use the Hong Kong island bus route data to 
analyze the impact of different parameters on the final result 
of BSDP optimization. CRO is also compared with SGA and a 
greedy method, SRM. Simulation results show that the proper 
selection of c has great impact on the final optimization result, 
namely, the total number of sensors needed, the area coverage 
percentage, as well as the data accuracy. Meanwhile CRO can 
outperform SGA and SRM with different values of c, in both 
average quality of solutions and the best generated solution. 

In the future we will conduct a systematic analysis on the 
variance of the different parameters and perform Student’s t- 
test on the variance. Then we will try to use metaheuristics 
other than CRO and SGA to see which metaheuristic is the 
most effective in solving BSDP. Moreover, the paper can be 
further extended in several ways. One possible extension is 
to have different number of sensors installed on different 
routes to reduce the sensing interval and increase the accuracy. 































Another is to set different coverage thresholds for different 
grids (regions) to reflect the different time constants for 
pollutant level changes. We will also deploy the sensors on 
a real transportation system for real-world testing. 
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